WHilSS  919 


UNCLASSIFIED 


THE  HORLD  SEODETIC  SVSTEH  1984  EIWTH  GMVITNTIONm. 
HODEL<U>  DEFENSE  NAPPING  A6ENCV  AEROSPACE  CENTER  ST 
LOUIS  AFS  NO  H  L  NHITE  82  NAV  8S 


I 


(.  NLLA'y'.  i » ;  lU 

.Hi  ■  Y*  <  A  I 'orj  <  .U  ImV 


0 


REPORT  DOCUMENTATION  PAGE 


'd  RtPOHI  SiCUKIIV  CLAiilUCAT'ON 
UNCLASSIFIEO  _ 


'  Zi  itCUR.rv  CLAiSlHCATION  AjTHOKiiY 

If)  N/A _ 


)b  RtSVKiCriVt.  MARKlNljS 
NONE 


3  OISTRIBUIION/ availability  'jf  RtPURT 


DTIC 


6c  ADDRESS  (Cty,  Stale,  and  ZIP  Code) 

St.  Louis,  Missouri 
:  63118-3399 


8a  NAME  OF  funding /SPONSORING 

8b  OfFiCE  SYMBOL 

ORGANIZATION 

(If  applicable) 

DMA  Aerospace  Center 

SO 

8c.  ADDRESS  (Of)',  Srafe,  and  ZIP  Code) 

St.  Louis,  Missouri 

63113-3399 

11  TiTlE  ilndudo  Security  CtJiSihcationI 

Tfie  World  Gcoaetic  System  1984  Earth  Gravitational  Model 


1i  I'fPSONAL  A'jTHORiS) 
Haschal  L.  Unite 


13a  type  of  report 
FINAL 


130  TiVIE  COVERED 
FROM  TO 


14  DATE  OF  REPOFIT  (Yea •,  Month,  Day)  IS  fA(j£  COUNT 
86  APRIL  28  -  MAY  2  25 


16  supplementary  notation 

Prepared  for  presentation  at  and  publication  in  minutes  of  Fourth  International  Geodetic  Symposium  on 
Satellite  Positioning  at  Austin,  Texas  April  28  -  May  2,  1986. 


1  / _ COSATI  CODES _ I  18  SUBJECT  TERMS  (Continue  on  reverse  if  necejiary  and  identify  by  block  number) 

~  F'Pto  I  GROUP  i  SL8  CPOyp  Geodesy,'  Earth  Ofavitational  Models!  Satellite^s'tioning .  • 


19  AaSTRACF  (Continue  on  reverse  if  necessary  and  identify  by  block  numberj 

The  World  Geodetic  System  1984  (WGS  84)  Earth  Gravitational  Model  (EGH)  consists  of  a  set  of 
normalized  geopotential  coefficients  complete  through  degree  (n)  and  order  (m)  180.  The  first  part 
of  the  EGM,  through  degree  and  order  41,  was  developed  as  a  weighted  least  squares  combination 
solution  fren  mean  freo-air  gravity  anomalies;  geold  undulations  developed  from  satellite  rada'' 
altimetry;  laser,  Doppler  and  NAVSTAR  Clobal  Positioning  System  (GPS)  satellite  tracking  data  and 
"lumped  coefficient"  data.  Procedures  used  in  the  EGM  development,  testing  and  evaluation  are 
discussed  with  particular  e-phasis  on  orbital  analysis  results  as  they  apply  to  Doppler  point 
posi tioning. 


OTIC  FILE  copy 


DO  FORM  1473,  84  VAR 


p  I  .4  MR  e'diOO”  m  jy  be  used  until  e»h  Jutted 
All  otn»  '  editions  ate  ohtole’e 


10  SOURCE  OF  r 

.'ND'NG  NUMBc'-/ 

PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO 

NO 

NO 

ACCESSION  NO 

ZO  DISTRIBUI  <JN  AVA'LABILUv  OF  ABSTRACT 

fB  LJNCl  ASS 'It'D  U'Nl'ViTfD  □  SAMI  aSRRT  □  DtiC  USERS 

Zl  abstract  security  CLASSIFICATION 

UNCLASSIFIED 

iya  NAM‘  O'  l•'-SPONS;|n,[  INU.V  UUA; 

HASCHAL  L.  WHITE 

ZZb  l  ELlRHQNE  (M</ucle  Are*  Co<*’) 

(  314)  263-4991 

iic  OFFICE  SYMBOL 

DMAAC/SDT 

_  Sf  (  i .K I T  V  r L ASSIF U  A  T ION  OF  'HIS  PAGE 

86  07  2 


NAME  AND  TITLE:  HASCHAL  L.  WHITE/GEODESIST 

ORGANIZATION:  Defense  Mapping  Agency 

Aerospace  Center 
3200  South  Second  Street 
St.  Louis,  Missouri  63118-3399 

TITLE  OF  ARTICLE:  The  World  Geodetic  System  1984 

Earth  Gravitational  Model 

WHERE  AND  WHEN  TO  Fourth  International  Geodetic  Symposium 

BE  PUBLISHED:  on  Satellite  Positioning;  The  University 

of  Texas  at  Austin;  Austin,  Texas; 

28  April  1986  -  2  May  1986 


SUSPENSE  DATE: 


15  APRIL  1986 


CERTIFICATIOM 


ORGANIZATION!  Defense  Napping  Agency 

Aerospace  Center 
Scientific  Data  Department 
3200  South  Second  Street 
St.  Louis,  Missouri  SSllS'SSDS 

TITLE  OF  ARTICLE!  The  World  Geodetic  System  1984 

Earth  Gravitational  Model 


WHERE  AND  WHEN  TO 
BE  PUBLISHED: 


REVIEWED: 


Fourth  International  Geodetic  Sysiposlum  on 
Satellite  Positioning;  The  University  of 
Texas  at  Austin;  Austin,  Texas;  28  April  1986  - 
2  May  1986 

/ Technical  Adequacy/Completeness 

Z  Does  not  contain  SAA  Data 

This  document  has  a  security  classification  of 

UNCLASSIFIED. 


I  have  reviewed  this 


abstract 


'^per  and  approve  it  for  presentation 


'Ajujt-iUhd 

*  ROBERT  H.  ^L  T 
Chief,  Scientific  Data 
Department 


hKfxaan  For 

NTIS  CRA&t 
DTIC  TAB 
Unannounced 
Justification 


□ 

□ 


By . . . 

Distibution/ 


Availability  Codes 


Avail  and/or 

Special 

THE  WORLD  GEODETIC  SYSTEM  1984 
EARTH  GRAVITATIONAL  MODEL 


HASCHAL  L.  WHITE 

DEFENSE  MAPPING  AGENCY  AEROSPACE  CENTER 
3200  SOUTH  SECOND  STREET 

.  ST.  LOUIS,  MISSOURI  63118-3399 

^  SUMMARY 

The  World  Geodetic  System  1984  (WGS  84)  Earth  Gravitational  Model 
(EGM)  consists  of  a  set  of  normalized  geopotential  coefficients  complete 
through  degree  (n)  and  order  (m)  180.  The  first  part  of  the  EGM,  through 
degree  and  order  41,  was  developed  as  a  weighted  least  squares 
combination  solution  from  mean  free-air  gravity  anomalies;  geoid 
undulations  derived  from  satellite  radar  altimetry;  laser,  Doppler 
and  NAVSTAR  Global  Positioning  System  (GPS)  satellite  tracking  data 
and  ^'fhmped  coefficient^^'  data.  Procedures  used  in  the  EGM  development, 
testing  and  evaluation  are  discussed  with  particular  emphasis  on  orbital 
analysis  results  as  they  apply  to  Doppler  point  positioning, 

1.0  INTRODUCTION 

The  procedure  used  in  the  development  of  the  WGS  84  EGM  through 
degree  and  order  41  was  to  form  normal  equations  for  each  of  the  various 
data  sets;  mean  free-air  gravity  anomalies,  mean  geoid'  undulations, 
Doppler  or  laser  tracking  data  from  each  satellite  used  and  "lumped 
coefficient"  data.  These  normal  equations  were  then  combined  one  at 
a  time  to  obtain  preliminary  or  intermediate  solutions.  All  of  these 
solutions  were  evaluated  by  comparing  differences  between  observed 
and  computed  mean  gravity  anomalies,  mean  geoid  undulations  and  the 
Doppler  residuals  from  selected  satellite  orbit  reductions.  The 
magnitude  of  and  changes  in  the  gravity  anomaly  degree  variances  computed 
from  these  intermediate  EGMs  were  also  carefully  monitored.  This 
procedure,  although  time  comsumlng,  ensured  prompt  identification  of 
any  problems  associated  with  a  particular  data  set.  Upon  completion 
of  the  degree  and  order  41  portion  of  the  WGS  84  EGM  using  least  squares 
techniques,  a  residual  1°  X  1°  (equiangular)  mean  gravity  anomaly  field 
was  developed  by  subtracting  the  contribution  of  this  EGM  from  the 
observed  1°  X  1°  field.  This  residual  mean  gravity  anomaly  field  was 
then  used  in  a  spherical  harmonic  analysis  to  extend  the  final  WGS  84 


EGM  to  degree  and  order  180.  This  expanded  WGS  84  £GM  produces 
significant  gravity  anomaly  and  geoid  undulation  differences  when 
compared  to  similar  values  from  the  degree  and  order  41  EGM  in  areas 
containing  short  wavelength  high  frequency  information.  On  the  other 
hand,  very  little  difference  occurs  in  areas  devoid  of  trenches  and/or 
ridges  or  in  areas  with  a  relatively  smooth  surface  gravity  field. 
For  moderate  to  high  altitude  satellite  orbit  computations,  the  degree 
and  order  41  portion  of  the  WGS  84  EOl  is  adequate.  Further  discussions 
in  this  paper  deal  primarily  with  the  tests,  evaluations  and  applications 
of  the  least  squares  derived  portion  of  the  WGS  84  EGM. 

2.0  DATA  USED  FOR  GEOPOTENTIAL  MODEL  DEVELOPMENT 

The  data  used  in  the  development  of  the  degree  and  order  41  portion 
of  the  WGS  84  EGM  consists  of  3°  X  3°  equal-area  mean  free-air  gravity 
anomalies,  3"  X  3°  equal-area  mean  altimetric  geoid  undulations,  data 
from  five  modern  and  two  historical  Doppler  tracked  satellites,  laser 
tracking  data  from  two  satellites,  GPS  satellite  tracking  data  (processed 
as  range  difference  data)  and  "lumped  geopotential  coefficient" 
information.  Each  of  these  data  types  and  the  processing  of  the  data 


into  normal  equations  are  discussed  separately. 

2.1  3”  X  3**  EQUAL-AREA  MEAN  FREE-AIR  GRAVITY  ANOMALY  DATA 

The  June  1984  Department  of  Defense  (DoD)  Gravity  Library 
equiangular  1°  X  1°  mean  free-air  gravity  anomaly  file,  balanced  with 
respect  to  the  WGS  84  Ellipsoidal  Gravity  Formula  and  with  terrain 
corrections  applied  (where  available),  along  with  their  error  estimates 
was  used  as  basic  input  to  compute  4584  3°  X  3°  approximately  equal-area 
mean  free-air  gravity  anomalies  and  their  corresponding  sigmas.  These 
gravity  anomalies  were  defined  in  terms  of  3°  latitude  bands  subdivided 
into  whole  degree  longitude  increments.  The  sigmas  used  as  input  to 
the  weighting  scheme  for  each  observed  3°  X  3°  mean  gravity  anomaly 
were  developed  directly  from  the  corresponding  sigmas  of  the  1°  X  1° 
data. 

The  mean  gravity  anomaly  normal  equations  were  formed  from  the 
observation  equations  which  included  all  nonzero  geopotential 
coefficients  below  degree  and  order  41  as  well  as  a  bias  parameter 
for  the  geopotential.  The  WGS  84  constants  for  the  earth's  rotation 
rate  (u),  the  ellipsoid  semimajor  axis  (a)  and  flattening  (f). 


and  the  product  of  the  earth's  mass  and  the  universal  gravitational 
constant  (GM)  %rare  used  in  these  equations.  Other  required  constants, 
such  as  those  necessary  to  define  normal  or  theoretical  gravity,  were 
coisputed  using  appropriate  WGS  84  constants.  The  weighting  scheme 
used  for  each  observed  3°  X  3°  equal-area  mean  gravity  anomaly  (1) 
is 

Wi  -  1.0/(oi*  +  o„*) 

where 

«  weight  of  the  i^^  3°X  3°  equal-area  mean  gravity 
anomaly 

■  sigma  for  the  1^^  3°X  3°  equal-area  mean  gravity 
anomaly  developed  from  the  1”  X  1°  sigmas 
Ogi  ■  an  a  priori  estimate  of  the  model  error  assumed 
for  each  3°  X  3°  equal-area  observation.  This 
represents  the  error  of  omission  resulting  from 
truncating  the  EGH  at  degree  41.  The  value 
assumed  was  ±10  milligals  (one  slgsia). 

2.2  3”  X  3”  EQUAL-AREA  ALTIMETRIC  MEAN  GEOID  UNDULATION  DATA 

A  1°  X  1°  equiangular  mean  geoid  undulation  data  file  developed 
from  SEASAT  radar  altimeter  measurements  was  merged  into  3°  -x  3°  approxi¬ 
mately  equal-area  mean  geoid  undulations.  Since  this  data  file  was 
latitude  bounded  as  well  as  restricted  to  oceanic  coverage,  it  consisted 
of  2918  3°  X  3°  mean  geoid  undulations.  An  additional  requirement 
when  forming  a  3°  X  3°  mean  value  was  that  at  least  two  thirds  of  each 
3°  X  3"  mean  geoid  undulation  be  oceanic  in  ocean/land  interface  areas. 
Error  estimates  for  these  Interface  areas  were  significantly  larger 
than  those  for  the  "fully  observed"  broad  ocean  area  geoid  undulations. 
As  in  the  mean  gravity  anomaly  normal  equations,  all  nonzero  geopotential 
coefficients  below  degree  and  order  41  and  a  bias  parameter  for  the 
geopotential  were  included  as  parameters.  Similarly,  WGS  84  constants 
were  used  throughout  when  computing  the  observation  equations.  The 
model  error  assumed  for  weighting  purposes  for  the  3°  X  3°  equal-area 
geoid  undulations  was  ±1  meter  (one  sigma). 

2.3  DOPPLER  DATA 


Significant  improvements  in  the  Doppler  satellite  tracking  network 
were  made  in  1971  when  the  collection  of  continuous  count  Doppler  data 


was  begun.  This  led  to  the  categorization  of  Doppler  tracking  data 
collected  before  and  after  1971  as  historical  and  modem  Doppler  data, 
respectively.  Other  equiinnent  changes  such  as  the  installation  in 
1975  of  rubidium  oscillators  in  the  fixed  Doppler  receivers  represent 
an  additional  improvement  of  the  modem  data.  Doppler  data  from  seven 
satellites,  including  five  in  the  modem  era,  was  included  in  the  WGS 
84  EGM  development.  The  common  names  and  orbital  characteristics  of 
these  satellites  are  tabulated  in  Table  1.  Two  six-day  data  spans 
were  selected  and  processed  for  each  of  the  seven  Doppler  satellites 
except  HILAT.  Due  to  time  limitations,  only  one  six-day  span  was 
processed  for  this  satellite. 

TABU  1 


DOPPLER  SATELLITE  ORBITAL  DATA 


SATELLITE 

NAME 

SEHIHAJOR 

AXIS  (KM) 

PERIGEE 

HEIGHT  (KM) 

INCLINATION 

(DEGREES) 

GEOS  3 

7216 

818 

115 

SEASAT 

7159 

778 

108 

NNS  68 

7457 

949 

90 

HILAT 

7179 

768 

82 

DB  14 

971 

63 

♦GEOS  1 

1113 

59 

♦BEACON  C 

7503 

941 

41 

♦HISTORICAL  SATELLITES 


These  data  spans  were  selected,  as  much  as  possible,  to  satisfy 
the  need  for  dense  Doppler  tracking,  good  balance  between  northern 
and  southern  hemisphere  tracking  stations  and  variation  in  the  arguments 
of  perigee.  These  basic  data  spans  were  very  carefully  edited  on  a 
points-within-a-pass  and  on  a  pass  basis.  After  the  data  editing  was 
completed,  normal  equations  were  formed  for  each  six-day  arc  consisting 
of  coordinates  of  all  data  contributing  tracking  stations  as  well  as 
for  all  geopotential  coefficients  through  degree  and  order  41.  Other 
parameters  included  in  the  normal  equations  were  those  for  satellite 
initial  conditions,  drag  multipliers,  radiation  pressure  and  solid 
earth  tidal  forces,  time  correction  parameters  (when  required)  and 
a  small  collection  of  analysis  parameters.  Bias  parameters 


were  mathematically  eliminated  as  part  o£  the  normal  equation  formation 
procedure.  The  Doppler  normal  equations  were  tested  individually  in 
E6M  tuning  solutions  and  in  various  combinations  as  part  of  an  extensive 
validation  effort  prior  to  their  incorporation  into  a  combined  set 
of  normal  equations  supporting  a  final  WGS  84  EGM  solution. 

2.4  LASER  DATA 

Two  laser  observed  satellites,  Starlette  and  LAGEOS,  were  selected 
for  WGS  84  EGM  exploitation  purposes.  The  laser  data  processing  and 
normal  equation  development  was  performed  by  the  Center  for  Space 

Research,  The  University  of  Texas  at  Austin,  under  a  Defense  Mapping 
Agency/Naval  Surface  Weapons  Center  contract.  The  Starlette  satellite 
with  its  SO  degree  inclination  and  805  km  perigee  height  supplemented 
the  Doppler  geodetic  satellites  for  the  determination  of  the 
geopotential.  The  data  span  selected  for  processing  was  the  93-day 
observational  phase  (August  through  October  1980)  of  the  MERIT*  Short 
Campaign.  During  this  period,  an  average  of  eight  passes  per  day  were 
collected  from  18  stations.  The  entire  93-day  data  span  was  treated 

as  a  single  dynamical  arc.  The  dynamical  model  used  included 
conventional  gravitational  and  solid  earth  tidal  forces^  a  60  term 
ocean  tide  force  described  by  Or.  Schwiderski  of  the  Naval  Surface 

Weapons  Center,  solar  radiation  pressure  and  the  Jacchia  1971  drag 
model.  Estimable  geodetic  parameters  included  the  same  geopotential 

coefficients  as  the  Doppler  medium  altitude  satellite  equations  together 
with  a  few  additional  resonant  terms.  Other  parameters  included  in 
the  laser  normal  equations  were  tracking  station  coordinates,  satellite 
initial  conditions  and  multipliers  for  the  radiation,  drag  and  solid 
earth  tidal  forces.  Bias  parameters  were  mathematically  eliminated 
in  the  process  of  forming  the  normal  equations. 

The  laser  tracking  data  available  for  the  LAGEOS  satellite  was 
sufficient  for  a  two  year  data  span  to  be  processed.  This  data  span, 
covering  1980  and  1981,  overlapped  the  MERIT  Short  Campaign  data  set 
chosen  for  Starlette.  The  LAGEOS  observations  represent  over  4000 
passes  of  data  collected  by  32  laser  tracking  stations.  As  with 
Starlette,  the  entire  data  span  was  treated  as  a  single  dynamical  arc. 

*MERIT  ■  Monitor  ^arth  Rotation  and  ^ntercompare  the  Techniques  of 
Observation  and  Analysis 


Estimable  geodetic  parameters  included  tracking  station  coordinates 
and  a  reduced  set  of  geopotential  coefficients.  However,  the  reduced 
geopotential  coefficient  set  was  selected  to  include  all  coefficients 
producing  orbital  perturbations  at  the  one  cm  level.  The  selection 
process  was  based  on  the  use  of  estimated  coefficients  an  order  of 
magnitude  larger  than  Raula's  rule  (10~^/n‘).  LAGEOS,  because  of  its 
high  altitude,  contributes  only  to  the  determination  of  the  lower  degree 
and  order  harmonic  coefficients. 

2.5  NAVSTAR  GLOBAL  POSITIONING  SYSTEM  (GPS)  DATA 

Four  continuous  weeks  of  simultaneous  tracking  data  from  five 
GPS  satellites  was  selected  as  the  GPS  data  set.  The  data  set  was 
characterized  by  the  absence  of  irregularities  in  the  atomic  clock 
time  histories  and  the  avoidance  of  eclipses  (the  entry  of  the  satellites 
into  the  earth's  shadow  causing  force  modelling  problems).  Separate 
sets  of  normal  equations  were  formed  for  each  week  of  data.  The  normal 
equations  included  50  geopotential  coefficients  as  parameters  and  allotted 
for  the  adjustment  of  the  universal  gravitational  constant,  a  systematic 
Z-axls  shift  and  a  scale  correction  for  the  tracking  network.  Other 
incidental  parameters  included  a  clock  frequency  and  aging  parameter 
for  each  satellite  and  station  in  addition  to  some  pass  and  station 
bias  parameters.  Due  to  their  altitude,  the  GPS  satellites  were  not 
expected  to  contribute  significantly  to  the  WGS  84  EGM.  However,  this 
data  was  included  because  it  could  possibly  result  in  a  slight 
improvement  of  the  WGS  84  EGM  for  GPS  orbit  applications. 

2.6  LUMPED  GEOPOTENTIAL  COEFFICIENT  DATA 

Lumped  coefficients  refer  to  certain  linear  combinations  of  zonal, 
low  degree  and  resonant  geopotential  coefficients  which  are  responsible 
for  rather  large  satellite  orbital  perturbations.  Analyses  of  fitted 
satellite  orbital  element  histories  yield  estimates  of  these  lumped 
coefficients  which  may  be  used  as  "observational"  data  in  determining 
the  applicable  geopotential  coefficients  themselves.  Numerous  papers 
have  been  published  on  this  subject  by  D.G.  King-Hele  and  C.A.  Wagner. 

An  extensive  literature  search  at  the  Naval  Surface  Weapons  Center 
produced  a  total  of  426  unique  observations  of  lumped  coefficients. 


Approximately  20  of  these  equations  were  deleted  Initially  on  the  basis 
of  the  author's  own  remarks.  Normal  equations  were  formed  with  the 
remaining  observational  values  and  standard  errors  taken  from  the 
literature.  The  parameters  of  the  nonsal  equations  generally  consisted 
of  all  relevant  geopotentlal  coefficients  through  the  41st  degree. 
An  exception  was  the  synchronous  satellite  data  which  was  truncated 
above  the  6th  degree.  Preliminary  validation  tests  were  not  completely 
successful  In  terms  of  duplicating  published  results  for  preliminary 
test  solutions.  However,  preliminary  test  solutions  did  produce  some 
residuals  ranging  from  five  to  15  times  the  standard  error.  This 
resulted  In  the  deletion  of  an  additional  20  observations.  Finally, 
the  remaining  observation  equations  were  combined  Into  two  sets  of 
normal  equations,  one  representing  the  observations  of  synchronous 
satellite  and  zonal  lumped  coefficients  and  another  representing  the 
higher  order  resonant  tesseral  geopotentlal  coefficients.  Unfortunately 
test  solutions  developed  with  the  higher  order  normal  equations  slightly 
degraded  rather  than  Improved  the  WGS  84  EGM.  As  a  result,  only 
synchronous  satellite  and  zonal  lumped  coefficient  observations  were 
Incorporated  Into  the  final  least  squares  portion  of  the  WGS  84  EGM. 
The  problem  with  the  higher  order  tesseral  harmonic  resonant  data  set 
has  not  been  Identified. 

3.0  WGS  84  EARTH  GRAVITATIONAL  MODEL  DEVELOPMENT  AND  EVALUATION 


The  WGS  84  EGM  through  degree  and  order  41  represents  the  solution 
of  a  set  of  normal  equations  formed  by  combining  the  Individual  normal 
equations  developed  from  the  previously  described  data  sets.  All  of 
the  Individual  normal  equations  were  combined  on  a  one  to  one  basis 
without  any  scaling.  The  Initial  step  was  to  combine  the  mean  gravity 
anomaly  and  geold  undulation  normal  equations  and  then  add  to  these 
combined  equations  the  satellite  normal  equations  one  at  a  time.  Each 
time  a  new  set  of  normal  equations  was  added  a  preliminary  EGM  solution 
was  obtained,  tested  and  evaluated.  Evaluation  was  accomplished  by 
computing  gravity  anomaly  degree  variances,  determining  the  differences 
between  computed  and  observed  3°  X  3°  equal-area  mean  gravity  anomalies, 
determining  the  differences  between  computed  and  observed  3°  X  3°  equal- 
area  mean  altlmetrlc  geold  undulations  and  by  analyzing  the  Doppler 


— ^  KT*  irw  »r«  v^  i^rT«r> 


residuals  obtained  In  orbit  reductions  of  a  six-day  GEOS  3,  a  four-day 
NNS  68  and  a  four-day  SEASAT  data  span.  All  of  the  data  spans  used 
in  the  orbit  reduction  analysis  were  distinctly  different  from  the 
data  spans  used  in  the  WGS  84  EGM  development.  Two-day  orbit  reductions 
were  also  accomplished  at  the  Naval  Surface  Weapons  Center. 

3.1  DEGREE  VARIANCES 

Gravity  anomaly  degree  variances  were  computed  for  each  Intermediate 
EGM  as  the  first  step  in  its  evaluation.  This  relatively  simple 
computation  serves  as  a  problem  indicator  when  larger  than  expected 
coefficient  magnitudes  occur.  Gravity  anomaly  degree  variances  are 
computed  by  the  equation 

n 


where 


m=0 

gravity  anomaly  degree  variance  in  mgal^  for  degree  n 


y  *  the  average  value  of  theoretical  gravity 

^n,m»  ®n,m  *  normalized  geopotential  coefficients  of  degree 
n  and  order  zs 

Gravity  anomaly  degree  variances  for  selected  geopotential  .models 
are  tabulated  in  Table  2.  The  first  of  these  models,  GAGH  (n=m>41), 
was  obtained  from  a  combination  of  the  3°  X  3°  equal-area  mean  gravity 
anomaly  and  the  3°  X  3°  equal-area  mean  geoid  undulation  normal 
equations.  Degree  variances  for  this  model  are  quite  similar  to  those 
obtained  for  the  WGS  84  (n*‘m*41)  EGM.  The  other  models  given  in  Table 
2  --  WGS  84  (n*m«36),  WGS  84  (n*m*30) ,  and  WGS  84  (n*m“24)  --  represent 
least  squares  test  solutions  of  the  WGS  84  EGM  combined  normal  equations 


with  the  parameter  set  limited  to  degree  and  order  36,  30  and  24, 
respectively.  The  degree  variances  for  these  test  solutions  agree 
quite  well  with  the  WGS  84  (n“m*41)  EGM  degree  variances. 

3.2  MEAN  GRAVITY  ANOMALY  COMPARISONS 

One  method  of  evaluating  an  earth  gravitational  model  is  to  compute 
the  mean  square  difference  between  mean  gravity  anomalies  developed  from 
geopotential  coefficients  (Agj,^  gravity  anomalies  developed 
from  observed  terrestrial  data  (Agt^-  ^he  terrestrial  field  used  for 
this  comparison  was  developed  from  observed  data  only  and  categorized 


TABLE  2 


GRAVITY  ANOMALY  DEGREE  VARIANCES 
UNITS  »  MILLIGALS* 


DEGREE 

GAGH* 

(n*in*41) 

WGS  84 
(n=m“41) 

WGS  84* 
(n»m*36) 

WGS  84* 
(n=m=30) 

WGS  84* 
(n*in“24) 

2 

7.6 

7.6 

7.6 

7.6 

7.6 

3 

33.6 

33.9 

33.9 

34.0 

34.0 

4 

19.2 

19.2 

19.2 

19.2 

19.2 

5 

20.8 

20.9 

21.0 

20.9 

20.9 

6 

17.6 

19.4 

19.4 

19.4 

19.4 

7 

20.6 

19.3 

19.4 

19.5 

19.5 

8 

10.2 

10.9 

10.9 

10.9 

9 

11.5 

11.5 

11.4 

10 

9.6 

11 

8.7 

12 

3.8 

13 

7.3 

14 

3.2 

3.2 

3.4 

15 

3.4 

3.4 

3.3 

16 

3.9 

3.9 

4.3 

17 

3.6 

3.6 

3.6 

4.3 

18 

3.0 

3.6 

3.5 

3.6 

3.7 

19 

3.2 

3.3 

3.4 

3.5 

3.6 

20 

2.4 

3.1 

3.2 

3.2 

3.5 

21 

2.4 

3.2 

3.3 

3.3 

3.5 

22 

3.4 

3.5 

3.6 

3.9 

4.3 

23 

2.5 

2.7 

2.8 

3.1 

Z.~8 

24 

2.3 

2.6 

2.6 

2.8 

2.2 

25 

2.9 

2.9 

2.9 

2.9 

26 

1.9 

2.4 

2.6 

2.9 

27 

2.0 

2.5 

28 

2.6 

3.4 

29 

2.6 

2.7 

30 

2.8 

3.0 

3.2 

31 

2.9 

2.9 

32 

4.1 

4.1 

33 

3.4 

3.9 

34 

4.6 

35 

4.2 

4.4 

4.4 

36 

2.9 

3.6 

2.8 

37 

2.8 

3.4 

38 

2.5 

2.8 

39 

3.1 

3.5 

40 

3.0 

3.6 

41 

2.5 

2.8 

*TEST  MODELS  FROM  LEAST  SQUARES  SOLUTIONS 


in  terms  of  percentage  of  observed  data  available.  For  example,  100 
percent  observed  requires  all  nine  of  the  1°  X  1°  equal-area  mean  gravity 
anomalies  within  a  3°  X  3°  equal-area  boundary  to  be  observed. 
Similarly,  six  out  of  nine  implies  67  percent  observed,  etc.  Equivalent 
logic  applies  to  the  percentage  of  observed  values  used  In  forming 
the  5°  X  5°  equal-area  mean  gravity  anomalies. 

WGS  84  EGM  through  degree  and  order  41  was  used  to  compute  3° 
X  3°  and  5°  X  5°  equal-area  mean  gravity  anomalies.  These  computed 
means  were  then  compared  to  their  observed  counterparts.  Comparisons 
were  also  made  at  different  truncation  levels  to  study  the  additional 
Informational  content  obtained  by  Increasing  the  degree  and  order  of 
the  EGM.  Other  models  evaluated  were  the  degree  and  order  24,  30  and 
36  EGMs  developed  as  test  solutions  of  the  WGS  84  EGM  data  set,  the 
WGS  84  (n*‘m*‘41)  EGM  without  the  geold  undulation  data,  and  an  EGM 
developed  from  a  mean  gravity  anomaly/mean  geold  undulation  only 
solution.  The  3°  X  3°  and  the  5“  X  5°  mean  gravity  anomaly  comparisons 
are  given  in  Tables  3  and  4,  respectively.  The  results  of  these 
comparisons  can  be  summarized  as  follows: 

a.  There  Is  no  appreciable  difference  between  the  results  obtained 
by  truncating  the  WGS  84  (n*m~41)  EGM  and  the  results  obtained  from 
the  truncated  solution  models.  This  indicates  that  the  aliasing  of 
higher  degree  and  order  Information  Into  lower  degree  and  order 
coefficients  Is  Insignificant  for  the  WGS  84  EGM. 

b.  The  difference  between  the  computed  and  observed  mean  gravity 
anomalies  decreases  as  the  degree  and  order  of  the  model  Increases 
Indicating  the  validity  of  the  higher  degree  and  order  coefficients 
of  the  WGS  84  EGM. 

c.  The  WGS  84  (n’'m’*41)  EGM  appears  to  be  almost  as  good  for 

representing  the  observed  mean  gravity  anomaly  field  as  the  model 
developed  from  the  mean  gravity  anomaly/mean  geold  undulation  data 
only. 

d.  The  Inclusion  of  mean  geold  undulation  data  In  the  WGS  84 
EGM  produced  a  model  that  agrees  better  with  the  observed  mean  gravity 
anomaly  field  than  the  model  developed  with  the  WGS  84  EGM  data  set 
less  the  geold  undulation  data. 


TABLE  3 


COMPARISON  OF  3“X  3®  EQUAL- AREA  MEAN  GRAVITY 
ANOMALIES  COMPUTED  FROM  EARTH  GRAVITATIONAL  MODELS 
WITH  THOSE  DERIVED  FROM  TERRESTRIAL  DATA 


'  '»■  '  1 

1 

EARTH 

GRAVITATIONAL 

MODEL 

• 

DEGREE 

OF 

TRUNCATION 

[<(Agt-Agh)*>]*5 

3311  OBS 
n  -  4007 

OBS 
n  -  3679 

100^  OBS 
n  -  3190 

WGS  84 

41 

±9.31 

±8.44 

±7.53 

36 

9.78 

8.99 

8.15 

30 

10.44 

9.73 

8.93 

24 

10.90 

10.26 

9.47 

WGS  84  EXPERIMENTAL  EARTH  GRAVITATIONAL  MODELS 


WGS  84 

36 

±9.85 

±9.03 

±8.15 

TRUNCATED 

MODELS* 

30 

10.56 

9.84 

8.95 

24 

11.03 

10.37 

9.50 

WGS  84  LESS 

41 

9.32 

8.52 

7".  79 

GEOID  UNDULATION 
DATA 

36 

9.82 

9.08 

8.33 

30 

10.51 

9.84 

9.10 

24 

10.97 

10.36 

9.60 

WGS  84  GRAVITY 

41 

8.80 

8.02 

7.34 

ANOMALY  AND 

GEOID  UNDULATION 

36 

9.35 

8.62 

7.97 

DATA  ONLY 

30 

10.12 

9.45 

8.79 

24 

10.67 

10.07 

9.38 

UNITS  -  MILLIGALS 

n  IS  THE  NUMBER  OF  SQUARES  INCLUDED  IN  THE  SAMPLE  OUT  OF  A  POSSIBLE 
4584  FOR  WORLDWIDE  COVERAGE. 

*LEAST  SQUARES  SOLUTIONS  OBTAINED  TO  THE  DEGREE  AND  ORDER  SPECIFIED 
USING  THE  WGS  84  EGM  DATA  SET. 


TABLE  4 


COMPARISON  OF  S^X  5®  EQUAL-AREA  MEAN  GRAVITY 
ANOMALIES  COMPUTED  FROM  EARTH  GRAVITATIONAL  MODELS 
WITH  THOSE  DERIVED  FROM  TERRESTRIAL  DATA 


EARTH 

GRAVITATIONAL 

MODEL 

DEGREE 

OF 

TRUNCATION 

f<(Agt-Agh)*>]*5 

40%  OBS 
n  -  1421 

80%  OBS 
n  *  1238 

100%  OBS 
n  -  1036 

WGS  84 

41 

±4.23 

±3.66 

36 

4.66 

4.14 

30 

6.67 

5.24 

4.63 

24 

7.27 

5.89 

5.20 

WGS  84  EXPERIMENTAL  EARTH  GRAVITATIONAL  MODELS 

WGS  84 

36 

±6.15 

±4.65 

±4.08 

TRUNCATED 

MODELS* 

30 

6.74 

5.33 

4.66 

24 

7.41 

5.98 

5.21 

WGS  84  LESS 

41 

5.83 

4.54 

4.13 

GEOID  UNDULATION 
DATA 

36 

6.20 

4.89 

4.47 

30 

6.76 

5.46 

4.94 

24 

7.34 

6.07 

5.44 

WGS  84  GRAVITY 

41 

5.11 

3.98 

3.52 

ANOMALY  AND 

GEOID  UNDULATION 

36 

5.53 

4.41 

3.94 

DATA  ONLY 

30 

6.21 

5.03 

4.48 

24 

6.92 

5.77 

5.13 

UNITS  •  MILLIGALS 

n  IS  THE  NUMBER  OF  SQUARES  INCLUDED  IN  THE  SAMPLE  OUT  OF  A  POSSIBLE 
1654  FOR  WORLDWIDE  COVERAGE. 

*LEAST  SQUARES  SOLUTIONS  OBTAINED  TO  THE  DEGREE  AND  ORDER  SPECIFIED 
USING  THE  WGS  84  EGM  DATA  SET. 


I  ^ 


3.3  GEOID  UMDULATION  COMPARISONS 

The  principles  involved  in  earth  gravitational  model  and  observed 
mean  gravity  anomaly  comparisons  can  be  extended  to  geold  undulations. 
However,  such  comparisons  are  limited  to  the  oceanic  geoid  determined 
from  satellite  radar  altimetry.  The  basic  data  set  is  1°  X  1° 
equiangular  mean  geoid  undulations  developed  into  3°  X  3°  and  5°  X 
5°  equal-area  mean  geoid  undulations  in  much  the  same  way  as  the 
equal-area  mean  gravity  anomalies  were  formed.  In  the  geoid  undulation 
comparisons  shown  in  Tables  5  and  6,  is  an  observed  mean  geoid 

undulation  developed  from  altimetry  and  is  a  mean  geoid  undulation 

computed  using  the  various  EGMs  being  evaluated.  The  results  of  these 
comparisons  can  be  summarized  as  follows: 

a.  As  in  the  mean  gravity  anomaly  comparisons  there  is  no 

significant  difference  between  the  results  obtained  by  truncating  the 
WGS  84  (n*m**41)  EGM  and  the  truncated  models  obtained  from  the  HGS 
84  EGM  data  set. 

b.  The  difference  between  the  observed  and  computed  mean  geoid 
undulations  decreases  as  Che  degree  and  order  of  the  model  used  to 
compute  the  geoid  undulations  increases. 

c.  The  models  developed  with  Che  geoid  undulation  data  included 

in  the  solution  produce  better  agreement  with  Che  observed  geoid 

undulations  than  the  WGS  84  EGM  developed  without  the  geoid  undulation 
data. 

d.  The  WGS  84  (n*m>41)  EGM  is  almost  as  good  for  representing 

geoid  undulations  as  the  model  developed  from  mean  gravity  anomaly/mean 
geoid  undulation  data  only. 

3.4  SATELLITE  ORBIT  ANALYSIS 

A  good  general  purpose  earth  gravitational  model  is  expected  to 
provide  a  means  for  computing  not  only  gravity  anomalies  and  geoid 
undulations  at  points  on  or  near  the  earth's  surface  but  also  to  serve 
as  adequate  gravitational  force  model  for  precise  satellite  orbit 

computations.  An  objective  in  Che  development  of  the  WGS  84  EGM  was 


TABLE  5 


COMPARISON  OF  3**  X  3°  EQUAL* AREA  GEOID  UNDULATIONS 
COMPUTED  FROM  EARTH  GRAVITATIONAL  MODELS  WITH  THOSE 
DERIVED  FROM  SEASAT  GEOID  UNDULATION  DATA 


EARTH 

GRAVITATIONAL 

MODEL 

DEGREE 

OF 

TRUNCATION 

WGS  84 

41 

36 

30 

24 

33Z  OBS 
n  -  3101 


t<(Nt-Nh)*>l>S 


67Z  OBS 
n  -  2918 


1007.  OBS 
n  -  2672 


±1.05 

1.16 

1.30 

1.46 


WGS  84  EXPERIMENTAL  EARTH  GRAVITATIONAL  MODELS 


WGS  84 

36 

±1.62 

±1.35 

±1.12 

TRUNCATED 

MODELS* 

30 

1.74 

1.49 

1.25 

24 

1.88 

1.67 

1.42 

WGS  84  LESS 

41 

1.99 

1.79 

1.59 

GEOID  UNDULATION 
DATA 

36 

2.03 

1.84 

1.63 

30 

2.12 

1.93 

1.71 

24 

2.22 

2.04 

1.82 

WGS  84  GRAVITY 

41 

1.47 

1.22 

0.96 

ANOMALY  AND 

GEOID  UNDULATION 

36 

1.55 

1.31 

1.07 

DATA  ONLY 

30 

1.67 

1.45 

1.21 

24 

1.83 

1.62 

1.38 

UNITS  -  METERS 

n  IS  THE  NUMBER  OF  SQUARES  INCLUDED  IN  THE  SAMPLE  OUT  OF  A 
POSSIBLE  4584  FOR  WORLDWIDE  COVERAGE. 

*LEAST  SQUARES  SOLUTIONS  OBTAINED  TO  THE  DEGREE  AND  ORDER  SPECIFIED 
USING  THE  WGS  84  EGM  DATA  SET. 


TABLE  6 


COMPARISON  OF  5°  X  5°  EQUAL- AREA  GEOID  UNDULATIONS 
COMPUTED  FROM  EARTH  GRAVITATIONAL  MODELS  WITH  THOSE 
DERIVED  FROM  SEASAT  GEOID  UNDULATION  DATA 


EARTH 

GRAVITATIONAL 

MODEL 

t<(Nt-Nh)*>]^ 

OEbKISE 

OF 

TRUNCATION 

40Z  OBS 
n  -  1108 

80%  OBS 
n  -  993 

100%  OBS 
n  -  887 

WGS  84 

41 

±2.34 

36 

2.35 

30 

2.41 

1.04 

24 

_ 

2.45 

1.46 

1.15 

WGS  84  EXPERIMENTAL  EARTH  GRAVITATIONAL  MODELS 


WGS  84 

36 

±2.35 

wmmm 

TRUNCATED 

MODELS* 

30 

2.41 

■B 

24 

2.50 

mgm 

1.12 

WGS  84  LESS 

41 

2.65 

1.67 

1.38 

GEOID  UNDULATION 
DATA 

36 

2.65 

1.69 

1.40 

30 

2.68 

1.74 

1.46 

24 

2.71 

1.80 

1.53 

WGS  84  GRAVITY 
ANOMALY  AND 

GEOID  UNDULATION 

41 

1 

36 

2.29 

2.31 

1.20 

1.23 

0.80 

0.85 

DATA  ONLY 

30 

2.35 

1.31 

0.94 

24 

2.40 

1.41 

1.07 

UNITS  •  METERS 

n  IS  THE  NUMBER  OF  SQUARES  INCLUDED  IN  THE  SAMPLE  OUT  OF  A 
POSSIBLE  1654  FOR  WORLDWIDE  COVERAGE. 

*LEAST  SQUARES  SOLUTIONS  OBTAINED  TO  THE  DEGREE  AND  ORDER  SPECIFIED 
USING  THE  WGS  84  EGM  DATA  SET. 


that  it  be  as  good  as  a  "tuned"  model  for  orbit  computations  where 
a  tuned  EGM  is  one  based  on  extensive  use  of  tracking  data  from  a  given 
satellite  in  the  EGM  development.  Within  this  context,  a  good  earth 
gravitational  model  for  orbit  computations  is  one  that  has  been 
extensively  tuned  by  including  satellite  tracking  information  from 
as  large  a  variety  of  satellites  as  possible.  However,  the  inclusion 
in  the  WGS  84  EGM  of  data  from  all  satellites  for  which  good  Doppler 
tracking  data  was  available  limits  orbit  analysis  as  an  independent 
method  for  evaluating  the  WGS  84  EGM.  A  second  consideration  is  that 
not  all  EGM  coefficients  produce  measurable  orbit  perturbations. 
Sensitivity  of  the  WGS  84  EGM  coefficients  through  degree  and  order 
41  at  or  above  the  0.5  meter  level  is  sho«m  in  Table  7  for  GEOS  3, 
in  Table  8  for  NNS  68  and  in  Table  9  for  all  ten  satellites  providing 
data  for  the  WGS  84  EGM  development.  The  blank  areas  in  Table  9  indicate 
those  coefficients  that  are  being  determined  primarily  from  mean  gravity 
anomaly  and  geoid  undulation  data.  Conversely,  these  coefficients 
are  not  being  evaluated  in  the  respective  orbit  analysis  tests  to  any 
degree  of  reliability. 

Data  spans  for  most  of  the  satellites  used  in  the  WGS  .84  EGM 
development  were  used  to  test  and  evaluate  the  WGS  84  EGM.  Orbit 
reduction  results  for  two  of  these  satellites  -*■  NNS  68  and  NOVA  -- 
are  presented  in  Table  10.  Data  spans  longer  than  the  period  of  the 
first  order  resonance  terms  were  selected  so  that  the  resonant 
coefficients  would  be  fully  tested.  The  tuned  EGM  for  both  the  NNS 
68  and  NOVA  satellites  (NWL  lOE-1)  was  used  for  the  initial  orbit 
reduction  and  data  editing  process.  These  data  sets  were  then  used 
in  orbit  reductions  with  the  WGS  84  (n*m«41)  EGM  without  any  further 
editing.  Since  the  purpose  of  these  tests  was  to  evaluate  the 
geopotential  coefficients,  other  constants  such  as  GM  and  tracking 
station  coordinates  are  common  to  both  the  tuned  model  and  the  WGS 
84  EGM  reductions.  One-day  drag  segments  were  used  for  NNS  68.  The 
NOVA  satellite  has  a  drag  compensation  system  that  in  effect  makes 
it  drag  free.  The  tabular  results  show  the  RWS  (the  square  root  of 
the  sum  of  the  squares  of  the  weighted  residuals  divided  by  the  sum 


WCS  84  CEOS  i  ORBITAL  PERTURBATIONS  >0.5  METER 


(^KEK  COEFFICIENTS  EXCLUDING  LUMPED  DATA  SETS 

2  lU  (ORBITAL  PERTURBATIONS  >0.5  METER) 


TABLE  10 


EVALUATION  OF  EARTH  GRAVITATIONAL  MODEL  COEFFICIENTS  IN  ORBIT  REDUCTION 
APPLICATIONS  (IDENTICAL  DATA  SETS,  STATION  COORDINATES  AND  GM) 


SATELLITE 

NAME 

EARTH 

GRAVITATIONAL 

MODEL 

RWS  (METERS) 

RAD 

TANG 

NNS  68 

NWT.  lOE-1* 

1.8 

2.8 

DAYS  205-208,  1982 

WGS  84 

1.5 

2.6 

NOVA 

NWL  lOE-1* 

1.7 

1.9 

DAYS  135-138,  1984 

WGS  84 

1.4 

1.9 

NOVA 

NWL  lOE-1* 

1.8 

2.0 

DAYS  141-144,  1984  ^ 

WGS  84 

1.6 

*  MODEL  TUNED  TO  THESE  RESPECTIVE  SATELLITES 


TABLE  11 

COMPARISON  OF  ORBIT  REDUCTION  RESULTS  BETWEEN  WGS  84  AND  THE  NSWC 
9Z-2  GEODETIC  SYSTEMS  FOR  THE  NOVA  SATELLITE 


GEODETIC 

SYSTEM 

_ om _ 

RWS  (METERS) 

DAYS 

YEAR 

RAD 

TANG 

NSWC  9Z-2 

135-136 

1984 

1.6 

1.6 

WGS  84 

1.5 

1.1 

NSWC  9Z-2 

137-138 

1984 

1.5 

1.5 

WGS  84 

1.2 

1.2 

NSWC  9Z-2 

141-142 

1984 

1.6 

1.7 

WGS  84 

1.6 

1.5 

NSWC  92-2 

143-144 

1984 

1.4 

1.4 

WGS  84 

1.3 

1.2 

of  Che  weights)  of  the  slant  range  (RAD)  component  (station-Co-satellite) 
and  the  tangential  (TANG)  or  intrack  component  of  the  Doppler  residuals 
resolved  at  TCA  (Time  of  Closest  Approach).  As  can  be  seen  from  the 
tabular  data,  the  WGS  84  EGM  produced  Doppler  residuals  as  low  as  and 
sometimes  lower  Chan  Che  tuned  models.  The  total  RWS  could  have  been 
reduced  in  all  cases  by  using  a  more  stringent  editing  criteria,  improved 
tracking  station  coordinates,  etc. 

Orbit  computations  were  also  made  using  NOVA  Doppler  Cracking 
data  for  two-day  data  spans  comparable  to  those  used  to  compute  precise 
ephemerides  for  Doppler  point  positioning.  In  these  computations, 
WGS  84  parameters  (EGM,  GM  and  station  coordinates)  were  used  for  the 
WGS  84  reductions.  The  NSWC  9Z-2  reductions  used  the  NWL  lOE-1  EGM 
and  GM,  and  NSWC  9Z-2  station  coordinates.  The  data  set  was  selected 
on  Che  basis  that  any  edited  pass  must  produce  large  residuals  in  both 
Che  WGS  84  and  NSWC  9Z-2  reductions.  The  results  of  these  computations 
are  given  in  Table  11.  In  all  cases,  the  Doppler  residuals  are  either 
equal  to  or  smaller  for  WGS  84  than  for  the  NSWC  9Z-2  geodetic  system. 

As  in  Che  case  of  Che  longer  data  spans,  extensive  editing  was  not 

required  since  the  objective  of  these  tests  was  to  evaluate  the  two 
systems  in  terms  of  Che  relative  magnitude  of  Che  Doppler  residuals 
rather  than  in  an  absolute  sense. 

3.5  WGS  84  EGM  CORRELATION  MATRIX 

The  correlation  matrix  for  Che  WGS  84  EGM  is  summarized  in  Table 
12.  This  Cable  indicates  those  coefficient  pairs  (C^j^oi  and  8^,01  for  a  given 
degree  and  order)  which  are  correlated  with  at  least  one  other  coefficient 
pair  at  a  level  >0.5  and  also  >0.7.  Comparison  of  Tables  9  and  12 
indicates  that  Chose  coefficient  pairs  that  produce  significant  perturba¬ 
tions  of  the  satellite  orbits  are  also  the  ones  with  correlation 

coefficients  >0.5.  This  relatively  high  correlation  is  as  expected  between 
the  satellite  sensitive  coefficients  because  of  the  similarity 
of  Che  period  of  the  orbital  perturbations  of  a  given  satellite  orbit 
for  coefficient  pairs  of  the  same  order.  This  higher  Chan  desired 
correlation  does  not  necessarily  present  a  problem  in  terms  of  some 
EGM  applications,  particularly  orbit  determination,  because  the  total 


or  "lumped"  effect  is  of  primary  importance  rather  than  the  singular 
effect  of  a  particular  coefficient  pair.  The  test  comparisons  and 
results  presented  for  both  mean  gravity  anomaly  and  orbit  computations 
indicate  that  the  WGS  84  EGM  is  not  significantly  affected  by  these 
higher  than  desired  correlations. 

4.0  CONCLUSION 

During  the  decade  since  the  development  of  WGS  72,  significant 
improvements  have  been  made  in  the  mean  gravity  anomaly  data  available 
for  gravity  modelling.  In  addition,  satellite  radar  altimeter  data 
for  determining  the  oceanic  geoid  has  also  become  available.  Other 
Improvements  such  as  accurate  Doppler  surveys,  improved  Doppler  tracking 
equipment  and  laser  tracking  technology  represent  significant  advances 
that  have  been  utilized  in  the  development  of  WGS  84  and  its  associated 
EGM.  The  test  results  and  comparisons  presented  here  demonstrate  that 
the  WGS  84  EGM  is  a  superior  EGM  for  DoD  applications.  Orbit  reduction 
tests  show  that  the  NNS  68  and  NOVA  orbits  used  for  Doppler  surveying 
can  be  improved  by  using  the  WGS  84  EGM.  This  improvement,  although 
small,  demonstrates  that  a  properly  developed  general  purpose  EGM  such 
as  WGS  84  is  as  good  or  better  than  the  "tuned"  EGM  previously  used 
by  DMA  for  NAVSAT  orbit  computations. 
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